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A theory of diffraction in the system consisting of the left-handed and the 

right-handed materials is proposed. The theory is based upon the Huygens's 

principle and the Kirchhoff 's integral and it is valid if the wavelength is smaller 

than any relevant length of the system. The theory is applied to the calcu- 
li ' 
00 

j — ' lation of the smearing of the foci of the Veselago lens due to the finite wave- 

o' 

' length. We show that the Veselago lens is a unique optical instrument for the 

^2 ' 3D imaging, but it is not a "superlens" as it has been claimed recently. 

In his seminal work Veselago has introduced the concept of the left-handed materials 

s 

(LHM's). In a simplest case the LHM's are materials with simultaneously negative electric 

a 

o 

o 



permittivity e and magnetic permeability \x in some frequency range. It is easy to show 
that in the LHM the vectors k, E, H form a left-handed set, while in the usual materials 



^ ■ (e > 0, fi > 0) they form a right-handed set. If imaginary parts of e and fi are small, 



the electromagnetic waves (EMW's) propagate in the LHM but they have some unusual 
properties. All these properties originate from the fact that in the isotropic LHM the 
Poynting vector S = E x H is anti-parallel to the wave vector k. 

Consider a propagation of the EMW from a point source located at the point z = —a 
through an infinite slab of the LHM with the thickness d and a usual right-handed material 
(RHM) at z < and z > d (Fig .[[]). It is obvious that S z > everywhere at z > —a 
because the energy propagates from its source. The directions of vector k for different rays 
are shown by arrows. They should be chosen in such a way that at both interfaces tangential 
components of vector k for incident, reflected and refracted waves are the same. Another 
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condition is that the component k z should be parallel to S z in the RHM and anti-parallel 
in the LHM. Then in the LHM k z is negative. It follows that the Snell's law for the RHM- 
LHM interfaces has an anomalous form: smi/ sinr = —n'/n, where i and r are the angles 
of incidence and refraction respectively, n' = -\/|e'||//|/eo/xo and n = \J e/i/e /io are positive 
refractive indexes for LHM and RHM respectively. The angles of reflection are equal to the 
angles of incidence. The refractive index in the LHM is often defined as negative 0, but we 
avoid this definition because in an infinite medium an EMW may propagate in any direction 
in both LHM and RHM and the only physical difference is that vector k is directed from 
the source of the wave in the RHM, while in the LHM it is directed toward the source. 

The device shown at Fig.|l|(b) is a unique optical lens proposed by Veselago. In this lens 
e = — e' and /i = — //, then n' = n and % = —r. It is easy to show that at n = n' the reflected 
wave is completely absent. Since all the rays going right from the source have % = — r, all of 
them have foci at points z = a and z = 2d — a as shown in Fig.|I|(b). 

All the ideas above have been put forward by Veselago about 35 years ago [I]]. Recently 
the method of fabricating of the LHM's on the basis of metallic photonic crystals has been 
found and the San Diego group has reported the first observation of the anomalous trans- 
mission and even the anomalous Snell's law Qj. Both observations have been interpreted 
as the result of negative e and \i. The speculations about the nature of negative e and // in 
photonic crystals are still controversial (compare 0-^]), but the very existence of the LHM 
seems to be demonstrated. Since the LHM's become reality it is time to develop a deeper 
understanding of their electrodynamic properties in order to use the advantages of these 
materials. 

One can see that the Veselago lens, shown at Fig.|l](b), is an absolute instrument because 
it images stigmatically a three-dimensional domain — d < z < and the optical length of 
any curve in the object space is equal to the optical length of its image 0. The only other 
absolute instrument we are aware of is the famous "fish-eye" of Maxwell . Note, that the 
definition of the absolute instrument assumes geometrical optics only. Since the LHM's have 
been already obtained we think that the Veselago lens can be extremely important device 
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for the 3D imaging. 

Pendry Q claims that the Veselago lens has a different unique property. Due to Pendry 
the resolution of this lens does not have a traditional wavelength limitation which follows 
from the uncertainty principle. Pendry has introduced a new term "superlenses" with the 



where the work of Pendry was criticized. 

In this paper we propose a general scalar theory of diffraction in the LHM which is based 
upon the Huygens's principle and the Kirchhoff's integral. As any diffraction approach 
our theory works under condition that the wavelength is much smaller than any relevant 
geometrical length in the problem. We apply this theory to the Veselago lens and calculate 
the smearing of the foci due to the finite wavelength. Thus, our result does not support the 
idea of "superlens" . The discrepancy between our result and previous ones is analyzed. 

The first problem is to find the Green function for the Helmholtz equation for the LHM 
which describes propagation of a spherical wave from the point source. It is easy to show 
that it has a form exp (—ikR)/R, where k = um/c and R is a distance from the source. At a 
small element of the sphere R = const the spherical wave can be considered as a plane wave 
which is characterized by the Poynting vector S and wave vector k both with the radial 
component only. Since S is directed along the external normal to the surface element, the 
wave vector k in the LHM is directed along the internal normal. It is easy to see that our 
Green function obeys these properties. 

Following the principles of the scalar theory of diffraction || the field u at the observation 
point P can be written in a form of a surface integral 



where R is the length of the vector from the point P to the surface element, df n is the 
projection of the surface element df on the plane perpendicular to the direction of the ray 
coming from the source to df, b\ is a constant for any LHM. To find bi one can consider a 
plane wave with the wave vector normal to the infinite plane of integration. Since this plane 



Veselago lens as a first representative of this class. Two comments appeared recently PJTH 




-ikR 



(1) 
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is fictional, the constant can be found from the condition that the Huygens's principle in the 
form Eq. ([]]) reproduces the same plane wave. Doing the calculations similar to Ref . || one 
gets bi = —k/2iri, so that for the LHM the constant bi has a different sign than the similar 
constant b r for the RHM. 

The Huygens's principle can be applied to any interface which has a curvature larger 
than the wavelength. It gives the correct direction of refracted waves but it does not give 
the amplitudes of both refracted and reflected waves. However, it can be successfully applied 
to the Veselago lens where reflected waves are absent. 

Note, that there are some other methods to describe the diffraction which may be also 
used if the source of the rays is unknown. They are described and compared in details in 
Jackson's textbook fLl] . O ne can show that all the methods give the same result at r = —i. 

Now we apply Eq.(|l|) to the Veselago lens. To find the field u inside the slab we shall 
integrate in Eq.([l]) over the plane z = 0. The field u in this plane is produced by a point 
source and has a form 

gifc^/ a 2 +x 2 +y 2 

u(x, y, 0) = =. (2) 

a/<2 2 + x 2 + y 2 

The field inside the slab can be found using Eq.([l|) with a constant b r i instead of bi because 
now we are integrating over the RHM-LHM interface rather than over the fictional surface 
in the LHM. In a similar way at the LHM-RHM interface one should use a constant bi r . 
Using the method described in Ref. || it is easy to show that b ri = b\ and b ir = b r . Thus 
one gets 



°T °T e iky/a 2 +xl+y'( e ~ik^/ z 2 +{xi-x) 2 +{yi-y) 2 

u(x,y,z) = ka / / — — — — 2 / o / x, , =^ dx idyu (3) 
J J a 2 + x{ + y{ ^z 2 + ( Xl - x) 2 + ( Vl - y) 2 
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where the additional factor a/ya 2 + x\ + y\ is the cosine of the angle between the ray, 
coming from the source to the point {xi, yi, 0} and the unit vector in z direction. One can 
see that the optical lengths for all rays (the sum of exponents in the integrand of Eq.(|3|)) 
from the point source to the focus, located at z — a, x = y = 0, are zero and the value of the 
field at the focal point u(0, 0, a) = ik, while the geometrical optics gives an infinite field in 



this point. To find u in the vicinity of the focus one should expand the integrand in Eq.([3D 
near the point (0, 0, a) assuming a; < a, 1/ <C a and \(\ <C a, where ( = z — a. One gets 



i sin (ky/ p 2 + ( 2 ) 



where p 2 = x 2 + y 2 . At p = one gets analytical expression 

1 — cos (k() + i sin (kQ 



J (kp\/l — s 2 ) sin (k(s)ds 
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Another analytical expression can be obtained at £ = 

i sin(A;p) 



u(p,0) 



(6) 



Figure ^| shows dimensionless function \u(p, C,)\ 2 /k 2 as given by Eq.(|4|). One can see that 
the smearing of the focus is anisotropic. The half-width in z direction is approximately one 
wavelength while in p direction it is approximately twice as less. At small x, y, ( the surfaces 
of a constant \u(x,yX)\ 2 are ellipsoids of revolution along z axis. 

Now we find the field u in the close vicinity of the second focus located at x = y = 0, 
z = 2d — a. The general expression for u at z > d differs from Eq.(|3|). One should apply 
the Huygens's principle to both interfaces located at z = and z = d. The later one is the 
LHM-RHM interface and the constant b r = k/2iTi should be used instead of hi. One gets an 
additional integral over the plane z = d so that expression for the field has a form 



u(x, y, z) = bib r ad 




dx\dy\ 




dxidyi 



a 2 + x\ + y\ 
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-iky/ d 2 + (xi-X2) 2 + (yi-y2)' 2 



-iky/ {z-d) 2 + (x2-x) 2 +(y2-y) 2 



(7) 



d 2 + (xi - x 2 ) 2 + (yi - y 2 ) 2 y/(z - d) 2 + (x 2 - x) 2 + (y 2 - y) 2 
To calculate these integrals in the vicinity of the second focus using inequalities kd ^> 1, 
ka ^> 1 one should introduce new variables {s,t} instead of variables {x2,y 2 } by relations 



x 2 = -{ l)Xi + s 

a 

,d 

y 2 = _( l)y l+t 



(8) 
(9) 



Equation (|j) has the following meaning at s = t = 0. For every point {xi,yi} at the 
first interface they give a point {£2,2/2} a ^ the secon d interface which is on the ray coming 
from {x\,yi} and passing through the first focus. Thus the new variables {s,t} describe 
deviation from the geometrical optics and they should be small. One can see from Eq.(^) 
that at s = t = optical lengths of all rays exiting from the point source at z = —a and 
coming to the second focus at z — 2d— a are equal to zero. Introducing the new variables and 
expanding the exponents in Eq.(^) one can get an expression for the field u in the vicinity 
of the second focus. For r\ = z — 2d + a, \rj\ <C a one gets 

u(x, y, rj) = -u(x, y, C)*\<;= v , (10) 

where the function u(x,y,() is given by Eq.(|3[). In other words, the smearing of the 
\u(x, y, r])\ 2 in the second focus is the same as the smearing in the first one. 

Equation (|To|) can be also obtained in a more physical way. One can calculate field u far 
from the foci expanding integrand in Eqs.(^|J^) near the geometrical rays as it was described 
above. In this case the result is exactly the same as in the geometrical optics. Namely, in 
the region < z < d 

ik(R-Ri) 

u{x,y,z) = R _ R ^ (11) 



where R = ay/1 + {x 2 + y 2 )/{a - z) 2 , R x = z^/l + [x 2 + y 2 )/{a - z) 2 . At z > d 



one 



should substitute into Eq.(|TT|) R = (z - d)y/l + (x 2 + y 2 )/(2d -a- z) 2 , R l = (d - 
a) yT+ (x 2 + y 2 ) / (2d — a — z) 2 . The expression in Eq.([lT]) becomes infinite in both foci 
as it should be in the framework of the geometrical optics. However, it can be used to 
calculate field u at the plane z = d 



g—iky/ (d—a) 2 +x 2 +y 2 

uix.y.d) — :■ (12) 

V ; y/(d - a) 2 + x 2 + y 2 V ^ 



Note that the minus sign in Eq.(|12|) results from the passing of the rays through the first 
focus. This gives extra phase tt. 

Now we can forget about the region z < d and apply the Huygens's principle to the 
z = d to find the field near the second focus. The field will be described by Eq. ([]]) with the 
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positive exponent and b r = k/2iri. Finally we get the result which is connected with Eq.([3D 
by Eq.©. 

Using Eq. ( |TTD one can calculate the flux of energy through any plane perpendicular to z 
axis for z > 0. One can show that it is independent on z and equal to 2n in our units. Note, 
that the flux of energy through the hemisphere around the point source at z = —a defined 
as f \u\ 2 df n is equal 2n, since u = exp {ikR)/R. 

Now we compare our results with the analytical calculations of Pendry [|J and Ziolkowski 
and Heyman fl2 |. Both papers claim that the Veselago lens in the ideal (lossless) regime is 



a "superlens" , which is able to provide a perfect focusing. In both papers the spherical wave 
outgoing from the source is represented as a superposition of plane waves, which are fictitious 
and do not correspond to the cylindrical symmetry of the problem. This superposition 
contains the "evanescent" waves (EW's), for which k 2 x + k 2 > u 2 /c 2 . One can easily show 
that the Poynting vector of each EW has a non zero components in the x-y plane but zero 
component in z direction. It follows from the second observation that the contribution of 
the EW's to the intensity near the foci of the Veselago lens should be exponentially small 
if ouj/c ^> 1. Pendry has explained perfect focusing as a result of amplification of EW's by 
the LHM. 



One can see from Eq.(47a) of Ref. |12| that the amplitude of a single EW increases 
exponentially in the LHM with increasing distance z from the source. Since the LHM is a 
passive medium, we think, that these EW-solutions should be omitted as nonphysical. The 
mathematical inconsistency of these solutions can be seen from the fact that the integral 
(Eq.(38) of Ref. [0), that describes the superposition of the plane waves, diverges at large 
k 2 + k 2 in the interval a < z < 2d — a at any value of x and y. Note, that the contribution 
of propagating waves into this integral is finite and it coincides with our result near the foci. 

The advantage of the diffraction theory is that it is a regular perturbation with respect 
to 1/ kd. We think that the EW's never appear in this theory because their contribution is 
of the order of exp (— 2kd). 



The computations, performed in Ref. [y|, do not show any focusing. We think that the 



main reason is that their a and d are of the order of the wavelength. Our calculations do 
not predict any focusing for such wavelength. 

Finally, we have proposed the theory of diffraction in a system, consisting of the LHM 
and the RHM and have applied this theory to the calculation of the smearing of the foci 
of the Veselago lens. This smearing is of the order of the wavelength so, from this point of 
view, the Veselago lens does not differ from any other lens. 

The work has been funded by the NSF grant DMR-0102964. 
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FIG. 1. Reflection and refraction of light outgoing from a point source at z = — a and passing 
through the slab of the LHM at < z < d. Refraction of light is described by the anomalous 
Snell's law. The arrows represent the direction of the wave vector. The reflected waves are shown 
by dashed lines near each interface only The slab is surrounded by the usual RHM. (a) n' > n. 
(b) The Veselago lens (n' = n). The reflected waves are absent, all rays pass through two foci. 




FIG. 2. Distribution of the dimensionless square modulus of the scalar field \u\ 2 /k 2 near the foci 
of the Veselago lens as a function of p and z as given by Eq.(||). Here \ = 2ir/k is the wavelength. 
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